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I. INTRODUCTION 



A consistent quantum theory of gravity should provide us with a clear picture of black- 
hole evaporation and black-hole entropy The main approaches to quantum gravity have up 
to now only led to partial successes []J , concerning mainly the statistical explanation for the 
entropy from microscopic degrees of freedom. 

In view of the fact that a generally agreed upon quantum theory of gravity does not yet 
exist, it is important to examine the quantization of particular models of physical interest 
by applying established quantization techniques. Among the most interesting models are 
those with spherical symmetry because they are sufficiently non-trivial and, at the same 
time, are simple enough to be tractable. One such model is that of Lemaitre, Tolman and 
Bondi (LTB), which describes a self-gravitating cloud of non-interacting dust. It was first 
introduced by Lemaitre (2J in an attempt to describe cosmology, where it has found inter- 
esting applications |3|]. It also describes gravitational collapse, but a full understanding of 
this process requires the application of quantum gravity. Given our limited understanding of 
quantum gravity at present, perhaps the approach that is most suited for addressing quan- 
tum collapse and questions concerning its final fate is canonical quantum general relativity. 
Although limited in its scope as a theory of quantum gravity, the canonical theory can 
meaningfully address some open questions regarding singularities in the quantum theory, 
such as the possibility of singularity avoidance, the quantum evolution of black holes, and 
the role of naked singularities in quantum gravity. 

In a recent paper we presented an exact canonical quantization of the LTB model 
Earlier work [si, @, 0, 0] had shown that the classical geometrodynamic constraints of this 
dust system, the Hamiltonian constraint and the momentum constraint, are both given 
in terms of a canonical chart consisting of the mass contained within spherical shells, the 
physical shell radius, the dust proper time and their conjugate momenta. The momentum 
conjugate to the mass function may be eliminated in the Hamiltonian constraint using the 
momentum constraint, and this procedure was shown to lead to a new and simpler constraint 
which is able to take the place of the original Hamiltonian constraint. Dirac's constraint 
quantization then yields the Wheeler-DeWitt equation for the wave functional describing 
the quantum collapse. 

In Q this approach was sharpened and particular attention was paid to the regulariza- 
tion of the Wheeler-DeWitt equation and diffeomorphism invariance. Regularization was 
performed on a one-dimensional lattice following a suggestion in [sj], and special care was 
taken to ensure that the momentum constraint is fulfilled in the continuum limit. As a 
consequence, the lattice wave functional becomes described in terms of three equations, the 
Hamilton-Jacobi equation and two additional constraints which make it possible to obtain 
exact solutions for a particular factor ordering. In this paper we use these exact solutions 
to examine Hawking evaporation during the collapse. The evaporation problem was also 
addressed by us in lOj , but only for marginally bound collapse in the WKB approximation. 
Here we consider the generic case, that is, the case including the non-marginal models for 
which the classical shells start with a non- vanishing velocity at infinity. 

Our paper is organized as follows. In Sec. II we review the major features of the quantum 
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LTB model. We present, in particular, the exact quantum states found in |4|. Sec. Ill is the 
main part of our paper. We show in detail how the spectrum of Hawking radiation together 
with appropriate greybody factors can be retrieved from these exact quantum states. This 
provides a bridge from quantum gravity (represented here by the LTB model) to semiclassical 
gravity and the Hawking effect. Sec. IV gives a summary of the obtained results and an 
outlook on potential future developments. 



II. THE CLASSICAL AND QUANTUM LTB MODEL 
A. The classical LTB model 

The LTB model describes a self- gravitating dust cloud. Its energy-momentum tensor 
reads T Mi , = e(r, p)u fl u u , where = u M (r, p) is the four-velocity vector of a dust particle 
with proper time r and labeled by p (p thus labels the various shells that together form the 
dust cloud). The line element for the LTB spacetime is given by 

ds i = _ dT 2 + (ffffi 2 dp i + R\ p )(d6 2 + sin 2 edcf) , (1) 
1 + 2E{p) 

and the Einstein equations lead to 

8nGe(T,p) = ^f^ and (d T R) 2 = | + 2E , (2) 

where F(p) = 2GM(p) is some non-negative function. We set c = 1 throughout. The case 
of collapse is described by d T R(r, p) < 0. 

There exists still the freedom to rescale the shell index p. This freedom can be fixed by 
demanding 

R(0,p)=p, (3) 

so that for r = the label coordinate p is equal to the curvature radius R. Now we can 
express the functions F(p) and E(p) in terms of the energy density e at r = 0. From (J2J) 
one gets 

F(p)=SnG [\(0,p)p 2 dp, (4) 
Jo 

E(p) = [d T R{r = 0, p)f --/' 6(0, p)p 2 dp . (5) 

P Jo 

The interpretation of these quantities is that F(p)/2G is the active gravitating mass inside 
of R(r,p), while E(p) is the total energy per unit mass of the shell labeled by p. The 
marginally bound models are defined by E(p) = 0. In the present paper we discuss the 
general case which includes the non-marginal case defined by E(p) ^ 0. 



The solution of (jzj) is given in [4[ and examined in great detail in [llj. The end state of 
collapse is either a black hole or a naked singularity, depending on the initial parameters 
governing the collapse. At a finite dust proper time r = r (p), the shell labeled by p reaches 
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a curvature radius R = 0, which results in a curvature singularity. Therefore r can only 
take values between — oo and r (p). 

The canonical formalism begins with the general ansatz for a spherically-symmetric line 
element, 

ds 2 = -N 2 dt 2 + L 2 (dr - N r dtf + R 2 dQ 2 , (6) 

where N and N r are the lapse and shift function, respectively. The total action is the 
Einstein-Hilbert action together with an action describing the dust. The phase space consists 
of the three functions L(t,r) and R(t,r) in (jHJ), the dust proper time, r(t,r), and their 
canonical momenta Pl, Pr, and P T . In the following, a prime will denote a derivative with 
respect to r. All variables are functions of t and r. After a series of canonical transformations, 
performed in the spirit of Kuchaf's reduction of static spherical geometries [5], it becomes 
possible to describe the phase space by a different chart, consisting of the dust proper time, 
r(r), the physical radius, R(r), the mass density function, T = F'(r) and their canonical 
momenta. The collapse problem is reduced to two classical constraints, 

H=(P T 2 + FP 2 )-^*0, (7) 

H r = t'P t + R'P R - YP T ' w , (8) 

where T = 1 — F/R, Pr is a new canonical momentum for R, which follows via a canonical 
transformation from the original variables, and Pr is the momentum conjugate to T. This 
transformation also absorbs a boundary term which is present in the original chart 



B. Relation between dust proper time and Killing time 



Under some circumstances, the momentum Pr has the interpretation of a Killing time. 
In this subsection we explain under which circumstances this interpretation makes sense. In 
the new variables, the proper time is constrained by the expression 



t' = 2PWI + 2E ± — Vl + 2E-F . 



T 



Defining a = 1/ yl + 2E gives 



(9) 



(10) 



a aT 

If the mass density vanishes for all r greater than a given r&, and if E is constant for all r 
greater than r&, Eq. (TTOT) can be integrated. This yields 



ar = 2P r ± / dR 



2Pr±P|4 Z | Z + ln 

1 - a 2 j2 



1 - Vl - a 2 F 



In 



1 + Vl - a?T 
Vl - a 2 ? - VT 



(11) 
(12) 



VI - o?T + VT^a? 
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where the upper sign refers to inf ailing dust and the lower sign to outgoing dust. Here we 
have assumed that < a < 1 (corresponding to E > 0). The Einstein equations guarantee 
that the discriminant 1 — a 2 T is non-negative. (A similar analysis can be carried out for 
-1/2 < E < 8). 

We know from Birkhoff's theorem that the spacetime around a collapsing dust cloud is 
given by the Schwarzschild solution. In |5] it was shown for the Schwarzschild geometry 
that 2Pr is equal to the Killing time T. Thus f ill I) connects the dust proper time r with the 
Killing time 2P T at the boundary r^. For small T and E' the relationship may still be used, 
since then we have a small amount of dust propagating in a Schwarzschild background; only 
in this case the concept of Killing time makes sense. This is the situation to be discussed in 
Sec. III. In the limit a — > 1 (E — > 0) we obtain 



r = T + 2\f¥ 



R 



(13) 



which is identical to the relation used in 10] for the marginal case. The plus sign after the 
T on the right-hand side has been chosen in order to describe a collapsing dust cloud. 

The transformation (TlT!) can also be interpreted in the following way. For a Schwarzschild 
spacetime it gives a relation between Killing time and the time used by families of freely 
falling observers. Each family is characterized by a fixed value of E. All observers within 
one family start at infinity with the speed v^, where E = f^ Q /(2(l — f^J). In the marginal 
case they thus start with zero initial speed. In this cas e ([TBI gives the relation between the 
Schwarzschild time and the Painleve-Gullstrand time [131 ]. 



C. The quantum LTB model 

Dirac's quantization procedure may be employed to turn the classical constraints in ([7]) 
and (E]) into quantum constraints which act on wave functionals. The Hamiltonian constraint 
(I7j) then leads to the Wheeler-DeWitt equation, and the momentum constraint (IE]) leads to 
a quantum equation which imposes diffeomorphism invariance on the wave functional. 

The resulting quantum constraint equations are of a purely formal nature. They must 
therefore be regularized before solutions can be obtained. Following a suggestion in ^ we 
have considered in Q] a one-dimensional lattice given by a discrete set of points rj separated 
by a distance a. In order that the momentum constraint is fulfilled in the continuum limit, 
it is important to start with a corresponding ansatz for the wave functional before putting 
it on the lattice. We therefore made the ansatz 

^ ]r(r), R(r), T(r)] = U drT(r)W(r(r), R(r), r(r))J , (14) 

where U : R — > C is at this stage some arbitrary (differentiable) function. The ansatz has to 
be compatible with the lattice, which means that it has to factorize into different functions 
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for each lattice point. So we have to make the choice U = exp, which gives 



^[r{r),R(r),T(r)) 

= lim TT exp^WiMn), R( ri ), F( n ))) 

i 

= lim TT %(r(n), R{n), r( n ), F( n )) , 



(15) 
(16) 

(17) 



(18) 



where 

i 

F(r,) = ^ar,. 

3=0 

The lattice version of the Wheeler-DeWitt equation was then found in 
where Tj = 1 — Fj/Rj > (outside the horizon). Making the ansatz ( 1171) together with 



to read 



*i = 



(19) 



the redefinition W = iW/2, we found in [4| that in order for the resulting equation to be 
independent of the choice of a (and thus also in the limit a — > 0) a set of three equations must 
be satisfied. If one finds solutions to all three equations, one can do all other calculations 
on the lattice, since these solutions have a well defined continuum limit and satisfy the 
momentum constraint. In ji| such solutions were found for a particular factor ordering in 
the Wheeler-DeWitt equation. (The solutions found in [9] do not fulfill the momentum 
constraints in the continuum limit.) 

The general solutions within the ansatz (fl4|) . with U = exp and satisfying the constraints, 
were found to be of the form 



^[t,R, T] =exp( J drb(F(r))T / Gh) 



x exp < =p 



2Gh 



dr T 



a(F(r))T± I dR- — 



(20) 



where we have introduced an implicit dependence a(F(r)) in order to get solutions fulfilling 
the momentum constraint, and b(F(r)) is an arbitrary function. They serve as the starting 
point for the discussion of the Hawking radiation in the next section. 



III. DERIVATION OF HAWKING RADIATION FOR THE GENERAL LTB 
MODEL 

A. Preliminaries 

The exact quantum states (120]) describe a generic situation. In order to study Hawking 
radiation in this framework, we need to introduce into the formalism the concept of a black 



hole and the analogue of the quantum fields used in the standard treatment. Following [10] , 
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we consider a Schwarzschild black hole plus some dust perturbation. The dust perturbation 



mimics the quantum fields in Hawking's derivation [14|. Formally, this perturbation is 
described by the variable r = F' in which the constant mass of the black hole drops out. 
In analogy to |Io| . the black- hole state factors out in (|20|) . and the remaining state will 
describe Hawking radiation. This remaining state then assumes the form of ( 1201) in which 
F can approximately be replaced by 2G times the black-hole mass. 



The central idea in [14|] was the use of a Bogoliubov transformation of field operators in 
the Schwarzschild spacetime. Here, the situation is different, since we have full quantum 
gravitational states at our disposal. We thus have to first identify those quantum states 
that correspond to the ingoing and outgoing modes, respectively, of the standard approach. 
Then we have to calculate their inner product. Since the description should refer to ob- 
servers at infinity, the inner product will be evaluated on hypersurfaces of constant Killing 
(Schwarzschild) time T, and not on hypersurfaces of constant dust time r (which corresponds 
to freely falling observers). 

We would therefore like to re-express our solutions in terms of T. For the expanding 
cloud we have the connection between dust time r and Killing time T given in ([12]) : 

ar = T- J R dR — ~ — . (21) 
The relevant exact solutions of the Wheeler-DeWitt equation read 



* ±=eXp (^/*'" rT 2^'"' 



CLT ± dR — 



(22) 



In principle one has to include additional normalization factors N + and iV~ for the wave 
functionals. Since it is not clear how to normalize the states properly in the full theory, we 
leave this point open here; it will not be crucial for our results derived below. 

In order to be as close as possible to the standard derivation of Hawking radiation, we 
want to obtain the scattered outgoing part of the 'ingoing wave functional'. We choose 
from fT22l) as in-going wave functional of positive frequency. This is justified as follows. 
Inserting the relation (121]) into the corresponding expression in (122]) . we obtain 



which has the standard form of a positive-frequency wave function. Note that this wave 
functional is independent of R. The lattice version of this state is obtained through the 
replacement of T = F' by a dimensionless quantity, u. This leads to 

+ _ Jj m TT ^-IbicrWiFilGh-ivUiTilGh 
i 

see f fTT]) . Our result will be a product of the results for each of the shells which together 
form the dust cloud. 
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For the out-going modes of negative frequency we have to take the state \I/ . Inserting 
( ]2Tj) into the corresponding state of (|22|) then gives 



cxp 



1 

Gh 



dr bT + 



2Gh 



dr T 



T — 2 



The corresponding lattice version reads 

2biOUJjFj iaujj 



lim I I 

(7^0 J- J- 



cxp 



J~ i 



(25) 



(26) 



In Sec. III.B we shall calculate the overlap of in- and out-going states in the near- horizon 
limit, while in Sec. III.C we give a general expression for the result away from the horizon. 



B. The Bogoliubov coefficient in the near-horizon limit. 

We shall now consider the analogue of the usual Bogolubov coefficient j3 uu i on each 
shell. We define it as the inner product between the in-going wave functional with positive 
frequency and the out-going wave functional with negative frequency (both frequencies are 
present in the Wheeler-DeWitt equation because of the P^-term there). In the standard 



approach one normalizes the field modes according to the Klein-Gordon inner product 15 
For solutions of the Klein-Gordon equation in 1 + 1 dimensions, 

<t> k (x,t) =u k (x)e- i ^ t , 

the spatial modes are normalized according to 

{u k ,u k ,} = — , 

2uj k 

from which one gets for the Bogoliubov coefficient 

/3 kk > = 2u k (u* k ,u k >) , 

where u k > denotes the 'new' set of field modes with respect to which the solution of the 
Klein-Gordon equation is expanded. We thus define in our case 

= ~Gn\ F dR y/todK&J ■ (27) 

Since we shall evaluate (3 separately for each shell, we skip the shell index % here and in 
the following. Note that au/Gh replaces the uj of the standard approach; it possesses the 
correct dimension of an inverse length, so that (3 itself is dimensionless, as it must be. The 
integration goes from the horizon {T = 0) to infinity. 

In order to specify the measure of integration, we recall again that the inner product 



should be evaluated on a T = constant hypersurface 10]. Performing the coordinate trans- 



formation (R, r) i — > (R, T) , we get with the help of ff2Tl) the result 



y/g^ = VT+2E^j; • (28) 



8 



We are now ready to evaluate the integral f|271) . Inserting the wave functional in (!!MI) and 
fl26|) into (J27D then gives 



Pu 



2uVl + 2Eexp (- 



2baF(uj + w') aTiuj + a/) 
i- 



Gh 



Gh 



dR^—^exp 



2iaw 
~~Gh 



(29) 



This is the exact expression, but it is too complicated to compute in the given form. However, 
in order to investigate black-hole radiation, a near-horizon approximation is sufficient 10 . 
Anyway, we shall derive in the next subsection an exact, albeit less tractable, result. In 
order to evaluate P, we introduce the variable 

r *-i 

F 

which ranges from zero (horizon) to infinity. The Bogoliubov coefficient, written in terms of 
s, then becomes 



P^< = 2wyl + 2i?exp 
2F 



2baF(uj + uo') .aT(w + w') 



Gh 



Gh 



X 



ds ^4— 4r exp 



o 



s 2 + 2s 



AiaujF 
Gh 



ds(l + s 



v /(l + 5 )2_ a 2 [(1 + 5) 2_ 1] 



s 2 + 25 



.(30) 



We have not specified the lower limit in the s-integral, since it will only contribute a phase 
to (3, which will not contribute to its absolute square. 

The near-horizon approximation is then described by an expansion with respect to s. 
This gives 



P, 



2Fau 

Gh 



Vl + 2Eexp (- 

POO 

/ dss- 1+2 ^ Gh exp 
Jo 



2baF(uj + to') aT(u + uj') 
i- 



Gh 

' auF3 + 10E ' 
1 Gh l + 2E S 



Gh 



x 



(31) 



In order to evaluate the integral in (13T]) . we insert a regularization factor exp(— ps), p > 0, 



which guarantees convergence at the upper limit. We use the formula [16 



dx x ^ e -^n> = T(v)(p 2 + g 2)-v/2 e -^arctan(g/p) 

and perform after integration the limit p — * 0. We then arrive at the following expression: 

2Fauu 



Pu 



IF aw , / 

_VTT2lexp (- 

(2iowF\ fawF3 + 10 E 
V Gh J \Gh 1 + 2E 



2baF(w + w') aT(w + w') 
i- 



Gh 

-2iau)F/Gh 



Gh 



x 



-iraujF/Gh 



(32) 
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For the absolute square of (5 we then get 



2 _ 47TF<TU(1 + 2E) e -*»F(uW)/Gh 

- ^ ^F/Gh _ i • ( 33 ) 

In order to calculate the particle-creation rate, we have to evaluate the expression J2k Pik/3*ki 
which here corresponds to the integral J °° du' \(3 UU) > | 2 (without further factors, in order to 
keep the result dimensionless) . We get 

H^m) = % 4 , wF/Cft _ 1 • (34) 

where iV out denotes the 'out'-particle-number operator. 

We recognize immediately that the spectrum contains a Planckian (thermal) factor and 
that the associated temperature is independent of E, as expected. In fact, the only in- 
dependent term is in the prefactor. 

Replacing auo by GAe, where Ae is the energy of a shell, we arrive at the final result 

,. , ~ ,. . 7rGAe(l + 2£) e~ 8bGMA ^ h 

<in|jV out |in> = eBKGMAe/h _ 1 > (35) 

where we have set F = 2GM. From here we read off the standard Hawking temperature for 
a Schwarzschild black hole: 

h 

k-aTu = , (36) 

which holds for each shell separately. 

We recognize that the Planckian spectrum is modified by greybody factors which explicitly 
depend on Ae. Greybody factors occur in the standard treatment by taking into account 



the back-scattering of field modes into the black hole; see, for example, [17|. Our greybody 
factors are different in origin compared to those coming from back-scattering, since we 
derived them from solutions of the Wheeler-DeWitt equation for gravity with quantum 
dust. 

At this stage, the number b has not been fixed further; it occurs in the full quantum 
gravitational state and depends on the imposed boundary conditions. The simplest choice 
for b is b = 0. In this case the integral over u' diverges. Dividing out the infinite constant 
arising from this integration, one arrives for b = at the result 

.. , ~ ,. , 8nGMAe(l + 2E) 1 

(m|iV out |m) = e s 7lG MAe/H_ 1 ■ ( 37 ) 



This corresponds to the result obtained in the WKB approximation, see 10| (after taking 
into account that the j3 there differs from the f3 here by a factor of 2a;). 



C. Exact expression for the Bogoliubov coefficient 

It is not difficult to obtain the exact expression for the Bogoliubov coefficient by consid- 
ering the series expansion in s about s = for arbitrary E of the integrand in fl30|) . The 
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measure can be expanded as 

\°- + ^ + ...^J2a n s«, (38) 



2(1 + s) 3 1 5 7s 



s 2 + 2s s 2 4 

n=— 1 



while the exponent becomes 



4iauF ( \ 1 . 9n 1 . 9 4n 9 

Ins + -(5 - 2a 2 )s + —(7 - 4a 2 - 2a 4 )s 2 + .... 



Gh \2 I 16 

def 2iacoF 

i ms + 

\ m=l 



Gh 



\ns + J2P m s m ) , (39) 



where we have returned to our notation a = 1/VT+2E. From pjj) we see that we must 
perform the integration 

n=-l 

where A = |(5 - 2a 2 ) = ±(3 + 10£)/(1 + 2E). This is best done by letting 
e 2 -^-^ m ^=l + ^f:/3 m s™ + i(^y P m P n s m s^jr im s™, 

m=2 ^ ' m,n=2 m=0 

(41) 

where all the coefficients 7 m can be evaluated in principle from the coefficients /3 m , and 
7o = 1, 71 = 0. We then find the integral 

°° /*oo 

= F <xn-ij m / dss- 1+n+m + 2i °" F ' Gh e 2i ™ F ^ s l GH 



n,m=0 



F V a n _ l7m f-i^-) / dss-^ F ' Gh e^ 
n j^ \ d X J Jo 



(42) 

X=2aujF/3 1 /Gh 



In this form, it can be evaluated as before. The first term (n = = m) is the near horizon 
approximation given earlier. The exact Bogoliubov coefficient can be written as the (infinite) 
sum 

2Fau , / 2baF(u + uj') aTiuo + uo 1 ) 

Pu** = -^-Vl + 2E exp I y — - i- 



Gh * V Gh Gh 

-2iaujF/Gh 

X 



f 2auF(3 i y— '™ F/Gh f 2iauF \ 

V Gh J ' V Gh J 

^Hr-(-¥)-(-¥--)S 

n+m=l \ / \ / \ / 



(43) 
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The simplest correction term to the near-horizon approximation is obtained for n = 1, 
m = 0: 

(-^{— G?r){-G?r) =-2£- 

Since 1 — 5/2/3i = —2/ (3 + 10-E), this leads to a cu-independent correction term in \/3\ 2 given 
by 

4 4 

< 



(3 + 10E) 2 ~ 9 ' 

In general, the ensuing correction terms in (143D depend on u; (and thus on the shell energy 
Ae) and therefore lead to a non-trivial modification of the greybody factors. In the limit of 
uj — > oo, the correction terms are independent of u>. 

Such a kind of greybody factors cannot be found in the standard derivation of Hawking 
radiation 141 ]. since there the geometric-optics approximation is strictly assumed. We can 
perform such a calculation because we have an exact quantum state at our disposal. 

We finally mention that no 'trans-Planckian problem' occurs here, since the whole cal- 
culation is based on genuine quantum gravitational states, and thus all Planck-scale effects 
are automatically included, at least within the limits given by our model. 



IV. SUMMARY AND OUTLOOK 

The programme for the canonical quantization of the LTB model was started in 
following methods earlier developed in [^j]. The classical Hamiltonian constraint and mo- 
mentum constraint were given in terms of a canonical chart consisting of the mass, radius, 
dust proper time, and their conjugate momenta. The momentum conjugate to the mass 
function was eliminated from the Hamiltonian constraint using the momentum constraint, 
and this procedure led to a new and simpler Hamiltonian constraint. The Wheeler-DeWitt 
equation for the wave functional describing quantum dust collapse then follows from the 
Dirac constraint quantization of the Hamiltonian constraint. 

To make progress towards finding solutions to the Wheeler-DeWitt equation, one must 
adopt a regularization scheme. The simplest choice is the delta- function regularization, 
which was used in jgj, but which turns out to be equivalent to the WKB approximation. 
The solutions for the marginal case in this approximation were obtained in j^] and were used 



to derive Hawking radiation in [lOj- In order to arrive at genuine quantum gravity effects 
one clearly needs to go beyond the WKB approximation. 

With this intent, a lattice regularization was proposed in @. A class of exact solutions to 
the Wheeler-DeWitt equation for the general non-marginal case were worked out in [4]. In 
order to construct such solutions, one assumes that the wave functional of the collapsing dust 
cloud can be written as a product of the wave functions of individual shells. We then look 
for solutions which are separable in r and R; it turns out that for the momentum constraint 
to be satisfied in the continuum limit, a specific factor-ordering must be chosen, which 
happens to depend on the energy function E. Interestingly, the separable exact solutions 
which then arise are a complete class of solutions, and there are no further non-separating 
solutions. Also, these solutions coincide with those in the WKB approximation. Thus in [4| 
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we found solutions for which the WKB form is exact - this is a consequence of the chosen 
factor ordering. This implies that the lattice regularization is in fact equivalent to the WKB 
approximation, with the additional constraint that the wave functional of the cloud be a 
product of shell wave functional. 

In the present paper we used these exact solutions to derive Hawking radiation and the 
accompanying greybody factors for the non-marginal case, thus generalizing our earlier work 
for the marginal case in 10j. Ideally, one would like to obtain quantum gravitational correc- 
tions to Hawking radiation - this, however, does not seem feasible within the framework of 
the presently used lattice regularization, since the WKB solution is also an exact solution, 
thereby ruling out any scope for a quantum gravity induced correction. (Interestingly, in 
18( | it has been suggested that such quantum gravitational corrections to the semiclassical 
Hawking radiation will be significant only above a certain minimum energy threshold). 

Another very important application of a model of quantum gravitational collapse is to 
investigate whether or not the singularity predicted for classical collapse is avoided in the 
quantum theory. In the present context this would require us to construct wave packets out 
of the exact WKB solutions that we have constructed here. Unfortunately, that does not 
seem possible for the non-marginal case since the factor-ordering we are constrained to choose 
depends explicitly on the energy function E (although one could perhaps construct wave 
packets with respect to b) . On the other hand, it may be possible to construct wave packets 
in the marginal case and we plan to investigate singularity avoidance in this case. Also, it 
would be of interest to investigate the LTB quantum collapse using the methods of loop 
quantum gravity, although in the loop approach one can at present strictly investigate only 
the cosmological singularity, and not the black hole singularity arising in an asymptotically 
flat spacetime. 

It also seems useful at this stage to look beyond the lattice regularization, and investigate 
other regularization schemes. The Wheeler-DeWitt equation for the LTB model is analo- 
gous to a two-dimensional quantum field theory. Hence it will be worth trying to adopt 
regularization methods used for two-dimensional quantum field theories to the present grav- 
itational context. We plan to take up such investigations for the quantum LTB model in 
the future. 
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